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Abstract 

We compute the DC conductance with two different methods, which both exploit 
the integrability of the theories under consideration. On one hand we determine the 
conductance through a defect by means of the thermodynamic Bethe ansatz and standard 
relativistic potential scattering theory based on a Landauer transport theory picture. On 
the other hand we propose a Kubo formula for a defect system and evaluate the current- 
current two-point correlation function it involves with the help of a form factor expansion. 
For a variety of defects in a fermionic system we find excellent agreement between the two 
different theoretical descriptions. 

PACS numbers: 73.50.Bk, 73.40.-c, 72.10.-d, ll.lOKk, 05.30.-d 



01alla@physik.fu-berlin.de 
Fring@physik.fu-berlin.de 



Contents 



1 Introduction 



2 Determining the defect scattering matrices 



2.1 Defect Yang-Baxter equations! 3 

2.2 Multiple defects! 5 

2.3 Free Fermion with defects 7 



p. 3.1 Tr ansmission and reflection amplitudes] 



2.3.2 



.3.3 



The energy operator defect V a {ip,ijj) 
Transparent defects, T >[) (t/j,ip) 



0, 



gipip 

2M 



.3.4 Energy insensitive defects, ^(ifj, ip) = gipj ^, T> 0± (ip,ip) = gi^{"f l ±7°)V 



.3.5 Luttinger liquid type V e (ifj,ip) = ipjgi + (feT^)^ 



.3.6 



.3.7 



The defects V^±{ijj,ip) 
The defects £> A± Q/>, -0) 



12 
12 
13 

~^7TT -')/2r| 14 

,#( 7 U ±-/ L )j2^ 14 



3 


Conductance from the Landauer formula 




3.1 Conductance through an impurity 




3.2 Defect TBA equations 




3.3 Thermodynamic quantities 






3.4 The high temperature regime 






3.5 Free Fermion with defects 




|3.5.1 Energy insensitive defects, T> u (ip,ip) = 0, T> IJ (tp , ip) , T>^ (tp , ip) , T> d± (x/),ip) 


|3.5.2 The energy operator defect T> a (ip , tp) = qtpip 



15 



20 
21 
21 



|3.5.3 Resonances versus unstable particles] 23 

|3.5.4 Multiple plateaux) 23 



4 Conductance from the Kubo formula 



24 



4.1 



4.2 



Conductance through an impurity 25 

Realization of the defect operator 28 

4.3 Defect matrix elements 29 

4.4 Free Fermion with defects 30 



4.4.1 



4.4.2 Conductance in the T 



Defect matrix elements! 30 



m 



4.4.3 Energy insensitive defects, V u (ip~ 



regime| 32 

o, ww^wws, WfWM 



4.4.4 The energy operator defect T>(ijj, ip) = gipip 



34 
35 



5 Conclusions 



37 



1 



1 Introduction 

Conductance (conductivity) measurements belong to the easiest and most direct experiments 
which can be carried out. They attract a lot of attention, due to the fact that in general 
they can be performed without perturbing very much the behaviour of the system, e.g. a 
rigid-lattice bulk metal, such that the uncertainty of experimental artefacts is reduced to 
a minimum. There exist various well-known theoretical descriptions, such as semi-classical 
transport theories (Landauer (l| and Boltzmann-Drude 0), dynamical linear-response theory 
II, H and also Green function linear-response theory 1| . To carry out the latter, in particular 
at finite temperature, is still poorly understood in generality ||, even in 1+1 space-time 
dimensions |7|. Since recent experimental progress allows conductance measurements also in 
1+1 space-time dimensions ||, one can on the theoretical side fully exploit the special features 
of low dimensionality. 

It is in particular very suggestive to exploit the full scope of non-perturbative techniques 
which have been developed in the context of integrable quantum field theories in 1+1 space- 
time dimensions, such as the thermodynamic Bethe ansatz (TBA) [0, O and the form factor 



bootstrap approach [11, 12]. Generalizing the Landauer transport picture a proposal for the 



conductance through a quantum wire with a defect (impurity) has been made in [p~3| , [LJ] 

oo 

G«(T) = £ lim | / de\pl(e,T,^\T?(9)\ 2 -pl(e,T,»l)\f?(e)\ 2 }, (1.1) 

1 — oo 

which we only modify to accommodate parity breaking, known to occur in integrable lattice 
models, see e.g. ]15| . This means in particular we allow the transmission amplitudes to be 
different for a particle of type i with charge qi passing with rapidity 9 through a defect of type 
a from the left T" (6) and right T-* (9). The density distribution function p\{9, T, //J, being 
a function the temperature T, and the potential at the left p\ and right ^ constriction of the 
wire, can be determined by means of the TBA. We have already restricted ( |1.1| ) to the abelian 
(diagonal) situation. It is clear that the effect resulting from the defect is most interesting 
when |T" (9) \ ^ 1, which requires the occurrence of simultaneous transmission and reflection 
(see Q, (1+7])). In this paper we will therefore be mainly interested in that situation. One 



may adapt fll.lj ) also to the case of pure reflection, which physically describes the influence of 
the constriction to the conducting process. From the previous statement it is clear that such 
boundary theories are only interesting in this physical context when they are non-abelian. 

The other prominent way of determining the conductance is a result from linear response 
theory, which yields an expression for the conductance in form of the Fourier transform of the 
current-current two-point correlation function. This Kubo formula has been adopted to the 
situation with a boundary 16]. As we mentioned, this will only capture effects coming from 
the constriction of the wire, we propose here a generalization to the analogous situation as 
described in (|l.l|) , i.e. when a defect is present 

oo 

G a (T) = - lim ^ J dt e*" (J(t)Z a J(0)) T ,m • ( L2 ) 



-oo 



Here the defect operator Z a enters in-between the two currents J within the temperature and 
mass m dependent correlation function. The Matsubara frequency is denoted by u. 
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The main purpose of this manuscript is to compare the two alternative descriptions (1.1 



and QljD for massive bulk theories with a defect which allows for simultaneous reflection and 
transmission. There exist various investigations, e.g., [17, O, 14, 18] for conformal (massless) 



theories with defect, which exploit the original folding idea of Wong and Affleck JlTJ. The 
idea is that a conformal field theory with a purely transmitting or reflecting defect can be 
mapped into a boundary theory, i.e. a theory living in half space, which has the advantage 
that the full restriction of modular invariance can be exploited in the construction of boundary 
states as pioneered by Cardy [19]. Since this folding idea relies on the vanishing of either the 
reflection or transmission, our considerations do in general not reduce to that set up, even in 
the conformal limit. As was already pointed out in [17], and as can be seen directly from (1.1) 
and ( |l.2[) , in that case the conductance is less interesting because it is either zero or perfect 
for abelian theories. 

In section 2 we outline the procedure of how the defect scattering matrices may be deter- 
mined, since they are needed as input in both approaches. In section 3 we newly formulate 
the defect TBA equations and use them to determine the density distribution functions. We 
evaluate numerically the Landauer formula ( |L1|) for various defects and provide some analyt- 
ical approximations in certain regimes. In section 4 we propose a Kubo formula ( |1.2| ) for a 
configuration in which an impurity is present and compute the current-current two-point cor- 
relation functions occurring in there by means of a form factor expansion. We find very good 
agreement between (|0|) and ( |l.2[) for the complex free Fermion theory with various types of 
defects. Our final conclusions and an outlook into open problems is provided in section 5. 



2 Determining the defect scattering matrices 

An essential input required in both non-perturbative methods which are exploited to compute 



the conductance (b_l) and (F2), that is the TBA and the form factor bootstrap approach, 
respectively, is the knowledge of the exact (defect) scattering matrix. It is one of the most 
intriguing facts of two dimensional quantum field theories that these matrices can be deter- 
mined exactly to all orders in perturbation theory. In the following section we will recall how 
much (little) of this approach can be carried over to the situation when defects are present 
and compute explicitly the transmission and reflection amplitudes for a variety of concrete 
defects. 

2.1 Defect Yang-Baxter equations 

A cornerstone in the context of integrable models in 1+1 space-time dimensions are the Yang- 
Baxter equations 123] . They can be derived most easily simply by exploiting the associativity 



of the so-called Zamolodchikov-Faddeev (ZF) algebra [21] and its extended version which 
includes an additional generator representing a boundary p2| , 23, ^4|] or a defect pq , |2q] . 
Indicating particle types by Latin and degrees of freedom of the impurity by Greek letters, 
the "braiding" (exchange) relations of annihilation operators Zi{6) of a particle of type i 
moving with rapidity 6 and defect operators Z a in the state a can be written as 

ZiiPdZjfo) = S%(6 1 -e 2 )Z k (6 2 )Z l (6 l ), (2.1) 
Zi{e x )Z]{e 2 ) = S%(6 1 -e 2 )zi(6 2 )Z l (6 l ) + 2ir5(8 1 -e 2 ), (2.2) 
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z0)z a = R^(e)z j (-e)z p + T?P(e)z p z j (e), (2.3) 
z a Zi(e) = Ri p a {-e)z p z 3 {-6) + f^{-e)z 3 {e)z p . (2.4) 

The bulk scattering matrix is indicated by S, and the left/right reflection and transmission 
amplitudes through the defect are denoted by R/R and T/T, respectively. We employed 
Einstein's sum convention, that is we assume sums over doubly occurring indices. We suppress 
the explicit mentioning of the dependence of Z a on the position in space and assume for the 
time being that it is included in a. For the treatment of a single defect this is not relevant, 
but it will become important when we consider multiple defects. The same relations hold 
when we replace the annihilation operators by the creation operators zj(9) with R/R, T/T 
and S replaced by their complex conjugates. The algebra ( |2.3] )-(|2~4]) can be used to derive 
various relations amongst the scattering amplitudes. Using extended ZF-algebra twice leads 
to the constraints 

6?6], (2.5) 

(2-6) 

0. (2.7) 

The same equations also hold after performing a parity transformation, that is for R <-> R 
and T <-> T in ( |2.6[) -( ^?7|) . From the associativity of the extended ZF-algebra one derives the 
equations ||, ||, |§ H II 

S(9 12 )[I ® i^(0i)]S(0 12 )[I ® R}{9 2 )} = [l® i^(0 2 )]S(0 12 )[I ® R^O^S^), (2.8) 
S(9 12 )[I ® RP^SihiW ® TJ{9 2 )} = flj(0 a ) ® Tg(0 2 ), (2.9) 
S(9 12 )[T£(9 2 ) ® T^(9 1 )} = [TP{9 X ) ®Tj(0 2 )]S(0 12 ), (2.10) 

where we employed the convention (A ® B)^ = A\B l - for the tensor product and abbreviated 

the rapidity sum 9\ 2 = 9\ + 9 2 and difference 9\ 2 = Q\ — 9 2 . Once again the same equations 
also hold for R <-> R and T <-> T. Starting with another initial asymptotic state one derives 

m 

R^M ® R}{9 2 ) = R}{9t) ® RP a (9 2 ), (2.11) 

[T^(9 2 )®l]S{h2)[R}{9i)®l]S{9 l2 ) = I%(0 2 ) ® Rtfa), (2.12) 

[l0T^(9 2 )]S(9 1 2W^R] i (9i)}S(9 12 ) = R%,{9\) ® fj{9 2 ), (2.13) 

[Ta{Qi) ® R]S(9x 2 )[Tp(9 2 ) ® I] = [l^f^9 2 )]S0 1 2)[l®T^{9 1 )]. (2.14) 

On the basis of the equations (|2^) -(^10|), it was shown in ]25t| , for the abelian case without 
defect degrees of freedom, that one can not have reflection and transmission simultaneously. 
In [^] this result was extended to the non-abelian parity breaking case and it was proven 
that for the simultaneous occurrence of reflection and transmission the scattering matrix has 
to be rapidity independent and of the form 

S(9)=¥a, (2.15) 



S%(9)S% n (-0) = 
R^(0)R-J(-0) + T/f(9)f ] k J(-9) = 

R^(e)T^{-e) + T^(9)R)}{-9) = 
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with P being a permutation operator and a a constant matrix. When assuming in addition 
that a is a diagonal matrix with the property (TijCji = 1, the free Fermion [pij = Gj{ = — 1), 
free Boson (cr^ = = 1) and also the Federbush model [27] and the generalized coupled 
Federbush models |2£| are solutions to ( 2.15| ). 



As a further set of consistency equations, which serve for the determination of the defect 
scattering matrix, we report the crossing relations, which are as usual less obvious to justify. In 
analogy to the relations which have to hold for the bulk scattering matrix Sij(9) = Sij{m—9) = 
0), (j is the anti-particle of j and * denotes the complex conjugation) we deduce from 



2^3) -(2.4) the crossing- hermiticity relations 

R°(6) = R°(-9)* = S j3 (29)Rf(in-9), (2.16) 
T?(0) = f*{-9)* = Tf(vn - 9) . (2.17) 

The first equalities follow when taking zj(9)* = Z{(9) and Z a = z\,. The latter relations in 
(2.17) simply result by considering the relations for S while letting one of the particles freeze, 



i.e., setting its rapidity to zero, and viewing it as a defect. Relations ( |2.16| ) are obtainable 
in a similar fashion as the interpretation put forward in |29|, p^j. Our equations ( |2.16| ) and 



( 2.17 ) disagree slightly from the crossing relations in ||29| , p?4| , |25| , pOf| , which is due to the fact 



that when parity is broken real analyticity is replaced by Hermitian analyticity [31|. Later 



on in our example, this will also be reflected in the representation of the free Fermion field 



( 2.35 ), being Dirac rather than Majorana. There is of course no consequence of this choice of 
conventions on the physics, since the ambiguity just exploits the fact that only the moduli of 
these amplitudes are observable. 

Similar as for the bulk scattering matrices an additional powerful constraint results from 
the singularity structure of the defect scattering amplitudes. In [^] it was shown that the 
defect does not admit any excited state once one demands a simultaneous occurrence of 
reflection and transmission. Supposing that the defect scattering matrices have a pole on the 
imaginary axis at i9q £ iR, the corresponding residues are therefore constraint as 

Res BS[0) = Res Rf(9) = Res T?(9) = Res ff{9) { < ° n f r ° a Q % ^\ . (2.18) 
e^ie 3 v ; 0->i0 o 3 e^ie 3 0^id o 3 \ > for 9 £ (0, tt) v ' 

The intervals (0, tt) are understood to be mod 2tt. Hence, there is no pole with positive residue 
in the physical sheet. 



2.2 Multiple defects 

Assuming that we have determined the defect scattering matrices R/R and T/T for a single 
defect, for instance by solving the consistency equations in the previous subsection, it is 
straightforward to use them in order to compute the related quantities for several defects. 
This type of situation is of interest since, unlike for a single defect, it leads to the occurrence 
of resonance phenomena and when the number of defects tends to infinity even to band 
structures. Let us therefore commence by exploiting the extended ZF-algebra Q2.3| )-( |2T4| ) for a 
double defect. For the reasons mentioned in the introduction we are interested in the situation 
when R/R and T/T are simultaneously non- vanishing, and in the light of the result ( p. 15 ), 
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we shall therefore focus on the diagonal case from now onwards. We compute 

Zi{9)Z a Z p = Rf(9)Z l (-9)Z a Z f3 + T^(e)Z a Z l3 Z l (9), (2.19) 

z a zpz { (9) = Bfi-e^aZpZi^ + Tfi-e^iWZaZp, (2.20) 

where we have now introduced overall transmission and reflection amplitudes corresponding 
to two defects 



T a P[e) = mwM R f {9)=Rm + Rf(o)TmTm (22i 

ff{9) = Tl^M , Rf( 9 ) = Rj(9) + ^ W^-gjfLgZ . (2.22) 



The term [1 - R P i {9)Rf(9)\ - 1 = £^° =1 (i?f (#)#? (#)) n results from the infinite number of 
reflections which we have in-between the two defects, well known from Fabry-Perot type 
devices of classical and quantum optics. For the case T = T, R = R the expressions ( [2.21 



and ( |2.22| ) coincide with the formulae proposed in [32]. When absorbing the space dependent 
phase factor into the defect matrices, the explicit example presented in [25] for the free fermion 
perturbed with the energy operator agree almost for T = T, R = R with the general formulae 
(2.21). They disagree in the sense that the equality of R^{9) and Rf l3 (9) does not hold for 
generic a, (3 as stated in [25]. 

It is now straightforward to extend the expressions to an arbitrary number of defects, say 
n, in a recursive manner 



rpai...a k ^Q^rpa k+ i...a n ^ 

l-Rf 1 - ak (9)R" h+1 - an {9) 



T, a (9) = K — i n l<k<n, (2.23) 

4 v ' 1 ijai...aw/i\ii a Hi- a »'«' 



fi° H1 - a " (9)T^- ak (8)f^- ak (9) 
l-Rf l - ak (9)R^ k+l - an (9) 



B»(0) = R^- ak (9)+ ' ^'l,^J;'^, n , > !<*<"■ (2-24) 



For convenience we encoded here the defect degrees of freedom into the vector a ={a\, • • • , a n }. 
Similar expressions also hold for ^(9) = f^- an {9) and Rf{9) = R^- a "(9). It is clear that 
from the knowledge of the single defect amplitudes we are now in the position to compute 
the corresponding quantities for multiple defects just by nesting successively the expressions 
(2.23) and (2.24) for increasing values of n into each other. Nonetheless, in general one does 
not succeed to provide simple analytical expressions for n-defect amplitudes and a different 
description is useful. 

Alternatively, we can define, in analogy to standard quantum mechanical methods (see 
e.g. |33|, [34|]), a transmission matrix which takes the particle from one side of the defect to the 



other. From the braiding relations (|2.3| ) and (2^), we obtain with the help of the unitarity 
relations ([2.6P and (| 



with 



(£///£52>) = (n *w>) * ^ 

M* (9) - ( T ? k(9rl -R? k (0)T? k (9)-i \ 
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This means alternatively to the recursive way ( |2.23| ) and ( p.24j) , we can also compute the 
multi-defect transmission and reflection amplitudes as 



\k=l J n \fc=l / i 2 \fc=l / n 

One may convince oneself that this formulation is indeed the same as ( |2.23| ) and ( |2.24| ). It 
has, however, the virtue that it allows for a more elegant computation of the band structures. 
In particular, it is most suitable for numerical computations, since it just involves matrix 
multiplications rather than recurrence operations. 

Let us now consider the case in which all the defects are of the same type a, equidistantly 
separated by an amount y and send n — > oo. First of all we have to include now explicitly the 
dependence of the defect on its position into the discussion. We assume 

" " if p ik y D \ 

HMi(x = ly) = H [Q y M l a {x = 0)] Q7 n ] , Q y = e _ lky , (2.28) 

i=i i=i ^ ' 

where k corresponds to the wavevector of the lattice. Taking then n — * oo this accommodates 
Bloch's theorem (e.g., ]33|] ) for the relativistic set-up. The simple requirement, that the 
product of transmission matrices limn-^tx, Y]!h=i Ma(x = ly) remains finite, leads now in the 
usual way to a restriction for the allowed energies, that is to band structures. To see when 
this is the case we can exploit the r.h.s. of the first equation in ( 2.2Sj ) and diagonalize the 



matrix QyM. % a {x = 0). Then it is clear that the limit n — ► oo only remains finite when the 
eigenvalues of this matrix are not real 



(2.29) 



The eigenvalues are computed to 



A l,2 



Xi(0) ± ^x?(0) 2 -f?(-d)/T«(-e), (2.30) 

x a {9) = [e^me)^ + e- ik y{ft{eyr l ] _ (231) 



In the parity invariant case the criterium (2.29) becomes simpler. From ( 2.30| ) and ( |2.31| ) 



follows in that case that the allowed energies in the infinite lattice have to respect 

Xi(9) = Re[e ifcy i; a (0)" 1 ] < 1, for T = f . (2.32) 

In other words particles are only allowed to travel in the system with rapidities for which 
the inequality (2.32) holds. In conclusion, this means the determination of the transmission 



amplitudes for a single defect is sufficient to determine multiple defects and the energy band 
structure. Let us illustrate the working of this general formulae with a concrete example. 

2.3 Free Fermion with defects 

The continuous version of the 1+1 dimensional free Fermion (Ising model) with a line of defect 



was first treated in [35]. Thereafter it has also been considered in |2^, g5[ and j30| from a 
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different point of view. In 



p?£| , |25[ the defect line has the form of the energy operator and 
in [3C] also a perturbation in form of a single Fermion has been considered. In this manuscript 
we want to enlarge the class of perturbations having in mind to obtain various different kinds 
of structural and physical behaviours. 

Let us consider the Lagrangian density for a complex free Fermion ip with I defects^ 



t-i 



£ = ipiij^df, - m)V> + S(x - x n )V a " {ifi,ifi) , 



(2.33) 



n=0 



where we describe the defect by the functions T> an (ip , ip) , which we assume to be linear in 
the Fermi fields tfi and tfi. In the following we will restrict ourselves mainly to the case of 
equidistantly distributed defects of the same type, i.e. x n = ny and V an (jfi, tfi) = Ufy,^) for 
n e {0,^-1}. 



2.3.1 Transmission and reflection amplitudes 

Unfortunately, it follows from the arguments outlined in section 2.1, that when one is seeking a 
situation with simultaneously occurring reflection and transmission the constraining equations 
for diagonal bulk scattering matrices reduce simply to unitarity and crossing. These equations 
are, however, not restrictive enough by themselves to fix R/R and T/T and therefore one has 
to resort to alternative arguments. For instance one may proceed in analogy to standard 
quantum mechanical potential scattering theory (see also |2^, [3(J) and construct the 
amplitudes by adequate matching conditions on the field. We consider now a single defect at 
the origin which suffices, since multiple defect amplitudes can be constructed from the single 
defect ones, according to the arguments of the previous section. We decompose the fields of 
the bulk theory as ifi(x) = @(x) ifi + (x) + 0(— x) ifi_(x), with 0(x) being the Heavyside step 
function, and substitute this ansatz into the equations of motion. This way we obtain the 
constraints 

&D(ip(x),ip(x))\ 



ij 1 (i' + (x) - ip_(x))\ x=0 



(2.34) 



x=0 



dip(x) 

Using here for the left (— ) and right (+) parts of ifi the Fourier decomposition of the free field 

dp] 



ifi Ax) 



a j (p)u j {p)e-^ x + aUp)v j (p)e' ip ^ x 



(2.35) 



with \ vrfc +p? = p°; and the Weyl spinors 




-0/2 
of/2 



and 




'Throughout the paper we use the following conventions: 
a; M = (a; ,^ 1 ), = (m cosh 8, m sinh 0), g 00 
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7 = 






-■=( 




. 7 



-6/2 

P e/2 



~e w = 1, 



(2.36) 



l t o 

■ 
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we can substitute them into the constraint (2.34). Treating the equations obtained in this 
manner componentwise, stripping off the integrals, we can bring them thereafter into the form 




f s (py r 3 (6)* 

T(9) RjiO) 
t. 




(2.37) 
(2.38) 



The creation and annihilation operators a{(9), a] (9) play in ( |2.1| ) and ( |2.2j ) the role of the ZF- 
algebra generators in view of the usual fermionic anti-commutation relations {aj(0ij, ajiOy)} = 
0, {ai(6i), a/jiOz)} = 2it5ij5{9i2). When including the defect operator in the equations ( |2.37|) 
and fl2.38| ), on the right/left for — /+-subscript, they acquire precisely the form of the extended 
ZF-algebra ( |2.3j )- (|2.4| ), such that one can read off directly the reflection and transmission 
amplitudes. One may convince oneself that the expressions found this way indeed satisfy the 

and respect ( |2.18| ). In 



consistency equations like crossing (|2.16| ), ( [2.17D , unitarity (|2.6D , (£ 
order to find the explicit expressions, we have to consider some concrete defects. Let us first 
concentrate on the energy perturbation. 

2.3.2 The energy operator defect V a (ip,ip) = gipip 

The defect T> a (^, tp) = g^pip has received already some amount of consideration, for the reason 
that it possesses a well studied |}6| discrete counterpart. Taking the continuum limit of these 
lattice models the defect term in there acquires the form of the energy operator e(x) = g^pip(x), 
with g being a coupling constant. According to ( 2.34} ), ( [2.37 ) and ( 2.38; ) we compute 



T«{9) 



Rf(e,y)=R?(e,-y) = Rf(e,-y) 

' {e) = Tf(o) = ff{p) 



sin B cosh I 



2iym sinh 8 



rpa 

3 



i sinh 9 
cos B sinh 9 



sin B 



sinh 9 + i sin B 



where we used a common and convenient parameterization in this context^] 

4 5 



sin!? 



7T 



4+<7 2 



< B < 0. 



(2.39) 
(2.40) 



(2.41) 



Note, that there is no explicit y-dependence in T/T and that ( |2.3S|)- ( 2.4Q|) sa tisfy the "unitar- 
ity" relations (|2.6| )- (|2.7| ) and the crossing-hermiticity relations ( |2.1€ )-( |2.17 ) when the defect 
is situated at the origin. The residues are constrained as in ( |2.18| ). The expressions R'j(9, B) 
and T°{6,B) coincide with the solutions found in |2j|, which, however, in general does not 
correspond to taking our particles simply to be self-conjugated, since we use Dirac Fermions. 
Having obtained these amplitudes, we can easily compute the corresponding quantities as- 
sociated to multiple defects by means of Q2.23 ), ( 2.24j ) or ( 2.27| ). The explicit formulae are 



f This is suggestive since many bulk theories admit such a relation between the bare and effective coupling. 
One may equate some combinations of R and T with some well known bulk scattering matrices. For instance, 
we identify the sinh-Gordon S-matrix Ssa(8, Bsa) — Tj(9,B-w/2)/Tj(—0,Bir/2), with the indicated relation 
amongst the effective defect coupling constants. 
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obvious and since they are quite cumbersome we will not report them here. Instead, we will 
depict them as functions of cosh# in figure 1 for various parameters in order to emphasize 
some of their characteristics. 




Figure 1: (a) Single defect with varying coupling constant. |T| 2 and |i?| 2 correspond to curves starting 
at and 1 of the same line type, respectively, (b) Double defect with varying distance y. (c) Double 
defect with varying effective coupling constant B. (d) Double defect = dotted line, eight defects = 
solid line. 

Part (a) of figure 1 confirms the unitarity relation (|2.6| ) where we used Rj(9,B) = 
Rj(—9,B) and T*(9,B) = Tj(—9,B). Part (b) and (c) show the typical resonances of a 
double defect, which become stretched out and pronounced with respect to the energy when 
the distance becomes smaller and the coupling constant increases, respectively. Part (d) ex- 
hibits a general feature which extends to an even number of higher multiple defects, say 2n, 
when keeping the distance y between the two most separated defects fixed: The resonances ac- 
cumulate at the position around the (2n— l)-th resonances of the double defect. For increasing 
n they become very dense in that region such that one may speak of energy bands. 

It is interesting to compare these bands with those obtained from the criterium ( 2.32| ), 
which translates in this case into 

sinh 9 (cos ky — cos B) < sin ky sin B, k = msinh.9 . (2.42) 
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Figure 1 part (d) shows that when taking 2n defects separated by a distance y/(2n — 1) one 
obtains for large n an energy spectrum which resembles a band structure. Analyzing instead 
the function X ?(0) in O) we obtain the same band structure from the criterium ( |2.32[ ) . The 
two computations show that the positions as well as the width of the bands in the two figures 
1(d) and 2 coincide quite well. Remarkably, even for the double defect the criterium ( |2.31| ) 
yields energy regions, see figure 2(b), which are in good agreement with the exact computation 
as presented in figure 1(d). 



A A A A A A A t 




Figure 2: Band structures according to the criterium ( |2.32| ). The non-shaded regions are forbidden, 
(a) eight defects with B = B\ = . . . = Bg = 1.1 equidistant by y = 0.25/7 (b) double defect with 
B = B 1 = B 2 = 1.1 distanced by y = 0.25. 

Very often we will not be able to perform certain computations analytically, but instead we 
can carry them out in the massless limit. The prescription for taking this limit was originally 
introduced in [37]. It consists of replacing in every rapidity dependent expression 9 by 9 ± a, 
where an additional auxiliary parameter a has been introduced. Thereafter one should take 
the limit a — > oo, m — > while keeping the quantity rh = m/2exp(o") finite. For instance, 
carrying out this prescription for the momentum yields p± = ±mexp(±#), such that one may 
view the model as splitted into its two chiral sectors and one can speak naturally of left (L) 
and right (R) movers. In this way the expressions ( |2.39D -( 2.40| ) become 

±ismBe ±2iyaThee 



R° 



,B) 



and 



L j,L/R\ u ^) - c anu T® L/R (9,B) = cos B. (2.43) 

Similarly we compute the expression involving two and four defects for later purposes 

cos 2 B 



rpa 1 a 2 i 



KB) 



n j,L/R\ 



rpaia2Ct-iCtA i 
1 j,L/R { 



,B) 



,B) 



,B) 



2iy a , me 



±ismBe- ziya i 



±ismBe 2iy ^ lheB ± 



1 + sin 2 B exp[^2im(y ai - y a2 )e e ] 
i sin £ cos 2 Se" 2 ^^ 



1 + sin 2 B exp[=F2im(y Q , 1 - y a2 )e 6 
ismB cos 2 Be~ 2iya ^ eB 



1 + sin B exp [^2iih(y ai - y a2 )e e ] 



T^ 2 ^ 3ai {9,B) 



rpCt\a 2 / 



R a \yj 



,B) 



(2.44) 
(2.45) 
(2.46) 
(2.47) 
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The remaining amplitudes can be obtained analogously. The expressions of physical quantities, 
e.g., the conductance, in the massless limit should not depend on the parameter m, such that 
the amplitudes ( 2.43| )-( 2.47| ) should in fact always appear in combination with other functions 
in order to make the prescription meaningful. 

Having discussed this type of defect in some detail we will now compute R/R and 1/1 
for various other defects in order to illustrate several types of physical behaviours. 



2.3.3 Transparent defects, P ^,^) = 0, V 13 ^),^) = g^ip 

The examples which can be handled most easily in later considerations are defects which 
behave physically as if they were transparent ones, i.e., as \T a \ = 1. Note that this does not 
necessarily mean the absence of the defect. For instance considering the defect T>" (^j) , ip) = 
gipj 1 ^, we compute with the method outlined above 



R 



Pi 



-B) 



.B) 



M(e,B) = o. 



-B) 



),B) 



JB 



(2.48) 
(2.49) 



for this defect. The coupling constant is parameterized as in ( 2.41 ). Evidently the "unitarity" 
( |2.6| )-( |2?7|) and the crossing relations ( 2.16| )-(2.17) are satisfied. Note that this is also an 
example for a defect which breaks parity invariance, i.e., the left and right transmission 
amplitudes are not identical. In the infinite lattice limit, i.e. when the number of defects 
tends to infinity, we find 



X PJ9) = cos(% T B) 



\Pje)£R,ve,B, 



which means that according to (|2.29| ) there are no forbidden energy regimes. 



(2.50) 



2.3.4 Energy insensitive defects, (tp , tp) = gifj^ip, T> s± (ip,ip) = gip(l l ±7 5 )'i/' 

In comparison with the transparent defects the next complication arises when the defect causes 

a phase shift independent of the energy of the incoming particle. For V 1 (ip , tp) = gipj^i/j we 

compute 



B,-y) 



R]( 



-B,y) = R]( 



-y) 



2iym sinh 6 



,B,y) = R]{6,-B, 

In this case we observe that parity is broken for the reflection amplitudes, i.e. R ^ R. The 



17(B) = 17(B) = 1](B) = 17(B) 



tan 5, (2.51) 
(2.52) 



relations (2.6)-( ^7[ ) and ( 2Tq )-( pT7 ) for y = are satisfied. For y = none of the amplitudes 
depend on the rapidities. In the infinite lattice limit we find 



X V$) = x](0) = cos ky cos B < 1 V0,B. 



(2.53) 



such that according to (|2.32| ) there are no forbidden energy regimes. 
For P <5± ( , 0,^) = g4>(^f l ± 7 5 )V' we compute 



R° 



,B,-y) = R 



-B,y)=R°±( 



1^(B) = T^i-B) 



T-. 



't(B) 



,B,y) = R S J ± ( 
1 5± (-B) = 



-B, —y) =±i tan — e 
B 

2i tan — . 
2 



2iym sinh 9 



,(2.54) 
(2.55) 
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These are examples in which parity is broken for the reflection as well as for the transmission 
amplitudes. Again the relations (|2^) -(l2?7|) and (|2^)-(|2~T?1) are satisfied when the defect 
is placed at the origin and, as for I? 7 , when y = none of the amplitudes depends on the 
rapidities. In this case we find in the infinite lattice limit 

^>)=^/ff) = Iif^ * \foO)tR,Ve,B, (2.56) 



such that according to (|2.32| ) there are no forbidden energy regions. 
2.3.5 Luttinger liquid type "D £ (t/>, ■(/>) = V>(<?1 +527°)^ 

When taking the conformal limit of a defect of the type T) 6 ^, ip) = ^(gi +527°)^ one obtains 



an impurity which played a role in the context of Luttinger liquids [38] when setting the 



bosonic number counting operator to zero, see e.g. |39|| . Besides T> a (ip , tp) this is also an 



example of a defect for which the potential is real. With ( 2.34 ), (|2.37|) and ( 2.38 ) we compute 
the related transmission and reflection amplitudes 



&m \ mm \ 4i(g 2 + gi cosh 9) e 2 ^ m ^ h 9 

Wn,9*-v) = R 3 (g. 9i ,92 ,y) = (4 - g2 _ g2) s . nh 9 _ 4z(gi - g2 cogh p , , (2.57) 

4i(si - g 2 cosh 9)e~ 2i y 171 sinhl 

Wm^-v) = Rj(0,9i,92,y) = (4 + g2 _ -^-^ _ ^ _ , (l'.w, 

rrsfn \ r^ein \ ( 4 + g| ~ gl ) sinh ,.-,r m 

T j (e, 9l , 92 ) = Tj (e, 91 , 92 ) = {4 + g 2_ g 2 )s . mhe _ 4l{gi + g2COshP , , (-'■*» 



\ x ( 4 + ffl -ffi)sinh6> 

(fl, ft , ft,) = I) (9, 91 ,92)= {4 + g2 _ gl) ginh g _ 4i(gi _ g2 cogh Q) ■ (2.60) 

As we expect, since lim 92 ^ ^ >e (V'> = '^ >a {i , -,'4')i we recover the related results also for the 
T/T's and ii/.R's in ( ^39|) - (^40|) . On the other hand, for 9l we obtain the defect 
T>^(tp,ip) = g^^ip f° r which the expressions simplify to 

_ • 2iymsinh6> : r> 

4(c?,S ; -y) = Rj(9,B,y) = — . -in,- ■ 5 u fl . ( 2 -61) 

■> J cos ±> sinn 9 + i sin is cosh # 

_ j p 2iym sinh • d 

i^((9, S, y) = RU9, B, -y) = -1— - — - — , (2.62) 



J ' cos -B sinh — i sin B cosh ' 



sinh 



Z?(0,B) = Tn9, B) = Tf(9, -B) = Tj(9 —B) = ^^- Q r^— 5 ttt-(2.63) 

J v ' ^ J v ' 3 ' cos B sum + isinB cosh 6> v y 



Where the effective coupling B is given by ( 2.41| ) with g — > 52 • The relations ( |2.6| )-(^ 



and (gT6D -(^l7|) may be verified once again for y = 0. In this case the infinite lattice limit 
leads to forbidden energy regimes, since according to ( |2.32| ), the rapidities have to respect the 
inequality 

(4 + g\ — g\) sinh# cos ky + sin ky{g\ ± g 2 cosh 9) < (4 + g\ — g\) sinh 6* for j, J, (2.64) 
which possess non trivial solutions. 
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2.3.6 The defects W± (i/j, if>) = gip(l ± 7 5 ) /2ip 
For this case we compute now 

^f9 —2iym s'mh9 

R]H0,B,y) - ^^-y)^ .^^^ (2.65) 

p±0 „2iym sinh 8 

r>, B) = ^ M |^W^M' lTl tM1 -, (g/2) si nh- w ■ (2 - 67) 

which is once again in ag reement with (|2~l|-(f27?l) and (^16|)-(^17|) for y = 0. In the infinite 
lattice limit we obtain also in this case forbidden energy regimes. The criterium ( |2.32| ) gives 

± cosky/2 < sinh 9 tan B/2 sin ky/2, (2.68) 

which has non-trivial solutions for the rapidities. 

2.3.7 The defects £> A ±(V>,V) = 9^{l° ±7 1 )/2i/> 
For this case we compute now 

—2iymsinh9 1. B_ 

R X + (9, B,y) = (9, B,-y)= , 6 g 2 , (2.69) 

J i sinh # — tan -j cosh 6/ 

— 2iym sinh j. _B 

^ (0, 5, y) = (0, B, —y) = . . , . f B , (2.70) 

J 2 sinn 9 + tan # cosh 9 

\, , , x,, , -x, , , -i, , , ± tan 5) sinh 9 

(9, B) = (9, -B) = lf± (-9, B) = T * (-9, -B) = , \ fl ^ B — ■ (2.71) 
J ^ « sinh 6> — tan ^ cosh # 

The crossing-hermiticity and unitarity relations hold for y = 0. 

In principle we could of course prolong this list of defects and construct their corresponding 
R's and T's. However, the main purpose of this exercise was to review how the transmission 
and reflection amplitudes for a defect may be obtained and also to give a brief account of some 
of their characteristic features. Important to note is that indeed all variations of possible parity 
breaking occur and one should keep therefore the discussion generic in that sense. Note that 
when the defect is real, namely T> a (ip,vp), T> £ (vp , ip) , parity invariance is preserved, which is 
a well known fact from quantum mechanics (see e.g. |54|). Complex potentials might look 
at first sight somewhat unphysical from the energy spectrum point of view. However, as is 
well-known for some bulk theories, such as for instance affine Toda field theories with purely 
complex coupling constants, one can still associate well defined quantum field theories to such 
Lagrangians and construct even classically soliton solutions with real energies and momenta 



ffl 



A classification scheme for possible defects which maintain integrability is highly desirable. 
It is interesting to note that in the conformal limit, as outlined before equation ( 2.43| ), some 



of the defects, namely T>^(t/),ip) and P A± (ip, i/j), become purely transmitting. Therefore, in 



contrast to first impression, the folding idea 17] could be employed. We have now enough 



examples at hand to use them in the following to determine the conductance in a multiparticle 
system, which we shall do in two alternative ways. 
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3 Conductance from the Landauer formula 
3.1 Conductance through an impurity 

The most intuitive way to compute the conductance is via Landauer transport theory Q. Let 
us consider a set up as depicted in figure 3, that is we place a defect in the middle of a rigid 
bulk wire, where the two halves might be at different temperatures. 




Figure 3: A conductance measurement. Part (a) represents the initial condition with no current 
flowing, i.e., 1=0 and part (b), I ^ 0. The defect is placed in the middle of the wire and the left and 
right half are assumed to be at temperatures T\ and T2, respectively. 

The direct current / through such a quantum wire can be computed simply by determining the 
difference between the static charge distributions at the right and left constriction of the wire, 
i.e. I = Q r — Qi- This is based on the assumption [13, 16], that Q(t) ~ (Q r — Qi)t ~ (p r — p h )t, 
where the ps are the corresponding density distribution functions. Placing an impurity in the 
middle of the wire, we have to quantify the overall balance of particles of type i and anti- 
particles 1 carrying opposite charges = —q^ at the end of the wire at different potentials. 
This information is of course encoded in the density distribution function p\ (9,T, //J. In the 
described set up half of the particles of one type are already at the same potential at one 
of the ends of the wire and the probability for them to reach the other is determined by the 
transmission and reflection amplitudes through the impurity. We assume that there is no 
effect coming from the constrictions of the wire, i.e. they are purely transmitting surfaces 
with T = T = 1. One could, however, also consider a situation in which those constrictions 
act as boundaries, namely purely reflecting surfaces. The situation could be described with 
the same transport theory picture, see e.g. 14, 41], but then the conductance can only be 
non- vanishing if the reflection amplitudes in the constrictions are non-diagonal in the particle 
degrees of freedom, such as for instance for sine-Gordon p9fl , that is in general affine Toda 
field theories with purely imaginary coupling constant or in the massless limit folded purely 
reflecting (transmitting) diagonal bulk theories. 

According to the Landauer transport theory the direct current (DC) along the wire is 
given by 



de 



P \{0,r,fi)\1? (0)| 2 - P mr,^\fr (6) 



(3.1) 
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d6 



pW t r,vZ)\R?(0)\'-pW t T,ri)\R?(e)\* , (3.2 



where we assume here T\ = T2. The relation (|3.2| ) is obtained from (3T) simply by making 
use of the fact that \R\ 2 + \T\ 2 = 1. Equation (|3.2[ ) has the virtue that it extracts explicitly 
the bulk contribution to the current which we refer to as 1b- There are some obvious limits, 
namely a transparent and an impenetrable defect 

lim I a = I B and lim I a = , (3.3) 

|T«|-»1 ' \T a \^0 



respectively. A short comment is needed on the validity of fl3.1| ). Apparently it suggests that 
when the parity between left and right scattering is broken, there is the possibility of a net 
current even when an external source is absent. In this picture we have of course not taken 
into account that charged particles moving through the defect will alter the potential, such 
that we did in fact not describe a perpetuum mobile. Thus the limitation of our analysis is 
that p\ — pj\ has to be much larger than the change in the potential induced by the moving 
particles. 

Finally we want to compute the conductance from the DC current, which by definition is 
obtained from 

G«(r) = ^.Gf(r) = ^. lim If (r, /(/4 - mD (3-4) 

and is of course a property of the material itself and a function of the temperature. In general 
the expressions in (|3.1|) tend to zero for vanishing chemical potential difference such that the 



limit in (3.4) is non-trivial. 

Thus from the knowledge of the transmission matrix and the density distribution function 
we can compute the conductance. Having already described how T" (9) can be determined, 
we will now explain how p^{9,r,p,j) can be evaluated by exploiting the integrability of the 
theory. 

3.2 Defect TBA equations 

The thermodynamic Bethe ansatz is a powerful tool which may be used to compute various 
thermodynamic properties of multi-particle systems which interact via a factorizing scattering 
matrix || and some chosen statistics. In addition, it allows to check the theory for consistency 
and to extract some distinct structural quantities such as the Virasoro central in the ultraviolet 
limit. The original bulk formulation has been accommodated to a situation which includes 
a purely transmitting defect [42] and a boundary [44|. In this section we want to propose a 



new formulation which is valid for a situation not treated before in this context, namely when 
reflection and transmission occur simultaneously. 

As usual we obtain the Bethe ansatz equation by dragging a particle along the world line 
of length L. We introduce for convenience the following shorthand notation for the product 
of various particle Zi(ff) and defect operators Z a 

Z ki',avM,a 2 ...k n ,a n := ^Mi(^Mi) • " " Z » kl (^ fcl ) Z "i • ■ ■ Z ^ kn ( G ^ k J Z a n ■ ■ ■ Z^d^). (3.5) 
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Then we compute the braiding of a particle operator of type i and the product of N further 
particles . . . Z Mjv with one defect Z a situated on the right of the particle Z^ k by using the 
ZF-algebra (O) and pj) 



Z t (6 t )Z^ 
Zi{9i) 



7 fi 1 ...fi N 

A k,a 



Zi(ei)z^ N F i( 



We abbreviated here 



pa 



N 



U s ^ e 



T?{-0i) ^ 
1 N 



till) 1 



N 



G1 



1=1 



1 l=k+l 
k N 



i=i 



n * 

=k+l 
V 

n SrtQ 



(3.6) 
(3.7) 



(3.8) 
(3.9) 



l=k+l 



Being on a circle of length L, we can make the usual assumption on the Bethe wavefunction, 
see e.g. M, which is captured in the requirement 



Zi(9)Z^ N = Z^ N Zi(e) exp(-iLm iS inh( 



(3.10) 



Using this monodromy property together with the braiding relations ( |3.6[) , (3/7) and the 
unitarity relation fl2.6|), we obtain 



N 

n 



\ SuiPu) 



N 



\i=i 



ALrrii sinh 9i 



T?{-Bi) e 



-iLrrii sinh£>j 



N 



X[Sii(Pu) 



(3.11) 



1=1 



Viewing the subscripts as entire spaces rather than components, equation ( |3.11| ) corresponds 
to the Bethe ansatz equation with simultaneously occurring transmission and reflection am- 
plitudes for the generic, that is also the non-diagonal, case. We restrict it here to the diagonal 
case and can therefore use the constraints ( 2.15| ), such that the bulk scattering matrix becomes 
rapidity independent and the relation (3.11) may be re-written as 



ALrrii sinh 9 



Su 



where 



DJfi) = 2 ± 2 



(3.12) 



(3.13) 

For consistency reasons it is instructive to consider the limit when the reflection amplitude 
tends to zero. In that case we can employ the relations ( |2.5|) -( ^?7|) and may take the square 
root in ( 3.13 ), such that we obtain from (|3.12| ) the two equations 



N 



N 



R,R^0 



ALrrii sinh 9 



ffWHSx, l = e 



iLrrii sinh 



X a (0)IlSu. (3.14) 



1=1 
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This means we recover the Bethe ansatz equations for a purely transmitting defect, which 
were originally proposed by Martins in [42]. The two signs in ( |3.13| ) capture the breaking of 



parity invariance in the limiting case, i.e. the two equations in ( |3. 14 ) correspond to taking 



the particle either clockwise or anti-clockwise around the world line as formulated for the 
parity breaking case for the first time in ||43||. We do not expect to recover from here the 



equations for a purely reflecting boundary which were suggested in [44], since the equations 
(|3.6| ) and (3/7) do not make sense in the limit T, f -> 0. For rjzli S% = 1, i.e. the free Boson 



and Fermion, we can exploit the fact that ( |3.12| ) with ( |3.13| ) look formally precisely like the 
Bethe ansatz equations for a purely transmitting defect. If we want to exploit this analogy we 
should of course be concerned about the question whether Dj a (6) is a meromorphic function. 
Assuming parity invariance, we may take the square root 

Df a [9) = Tf(d) ± R$(0) for T = f,R = R. (3.15) 

The matrix D^ a {9) has now the usual properties, namely it is unitarity in the sense that 
Dfj6)Df a (-d) = 1. It follows further from (gTp , (fyj) and (|2T7|) that the hermiticity 



relation Df a {9) = Df a {-9)* and the crossing relations Df a {9) = Dj a (iir - 9) and Df a (6) = 
Dj a (iTr — 0) hold for the free Fermion and Bosons, respectively. 

Let us now carry out the thermodynamic limit in the usual way, namely by increasing 
the particle number and the system size in such a way that their mutual ratio remains finite. 
The amount of defects is kept constant in this limit, such that there is no contribution to 
the TBA-equations from the defect in that situation, see also ^| where the same argument 
was employed. Hence this means that essentially we can employ the usual bulk TBA analysis 
when the considerations are carried out not per unit length. 

Let us therefore recall the main equations of the TBA analysis. For more details on the 
derivation see and in particular for the introduction of the chemical potential see [1C]. 



The main input into the entire analysis is the dynamical interaction, which enters via the 
logarithmic derivative of the scattering matrix (fij(9) = —idlnSij(9)/d9 and the assumption 
on the statistical interaction, which we take to be fermionic. As usual [^, [l(|, we take the 
logarithmic derivative of the Bethe ansatz equation ( |3.12| ) and relate the density of states 
Pi(0, r) for particles of type i as a function of the inverse temperature r = 1/T to the density 
of occupied states p\{Q,r) 

Pi (9,r) = T ^co S h9 + J2 j l<PiJ*Pim- (3-16) 

By (/ * 9) ($) := V(2tt) / dO 1 f(9—9')g{9) we denote as usual the convolution of two functions. 
The mutual ratio of the densities serves as the definition of the so-called pseudo-energies 
Si(9,r) 

Pi (e,r) l + e-^y { ' 

which have to be positive and real. At thermodynamic equilibrium one obtains then the 
TBA-equations, which read in these variables and in the presence of a chemical potential Pi 

rmi cosh 9 = Ei(9, r, Pi ) + + ^ . [tp^ * ln(l + e~ E ' )](6>) , (3.18) 
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where r = m/T, mi — ► mi/m, fa — ► fa/m, with m being the mass of the lightest particle in 
the model. It is important to note that fa is restricted to be smaller than 1. This follows 
immediately from ( |3.18| ) by recalling that £j > and that for r large Ei(6,r,fa) tends to 
infinity. As pointed out already in Q (here just with the small modification of a chemical 
potential), the comparison between (3.18) and (|3.16| ) leads to the useful relation 

Pi(0,r,Hi) = — [ + fa). (3.19) 



2tt 



dr 



The main task is therefore first to solve (3.1c: 
densities can be reconstructed. 



for the pseudo-energies from which then all 



3.3 Thermodynamic quantities 

Treating the equations ( 3.12 ) and ( 3.13j ) in the mentioned analogy we can also construct 
various thermodynamic quantities. It should be stressed that these quantities are computed 



per unit length. Similarly as the expression found in [42] for a purely transmitting defect the 
free energy is 



Fir) 



7TV 



La 



mi 



dO [cosh + m 1 <fi a (0)] ln[l + exp(— rm cosh ( 



(3.20) 



It is made up of two parts, one coming from the bulk and one including the data of the 
defect in form of <Pi a (0) = —id\nDi a {0)/d9. From equation ( 3.2C| ) we also see that when 
taking the mass scale to be large in comparison to the dominating scale in the defect, the 
latter contribution to the scaling function becomes negligible with regard to the bulk and vice 
versa. 



3.4 The high temperature regime 

Since the physical quantities require a solution of the TBA-equations, which up to now, due 
to their non-linear nature, can only be solved numerically, we have to resort in general to a 
numerical analysis to obtain the conductance for some concrete theories. However, there exist 
various approximations for different special situations, such as the high temperature regime. 
For large rapidities and small r, it is known ]g] (here we only need the small modification of 
the introduction of a chemical potential fa) that the density of states can be approximated 
by 

m, : i fl i 1 , dei(6,r,fa) 



Pi(8,r, fa 



-m- 



where e(0) = 6(0) 



47r 2irr 
0(— 9) is the step function, i.e. e{9) 



dO 



(3.21) 

1 for 9 > and e(0) = -1 for 



< 0. In equation (|3.17| ), we assume that in the large rapidity regime p[(0, r, fa) is dominated 



by (3.21) and in the small rapidity regime by the Fermi distribution function. Therefore 

p r i (6,r,fa)~^e(9)-^\n[l + eM-£i((>,r,fa))] • (3.22) 
~~ D, we approximate the direct current in the ultraviolet 



Using this expression in equation (| 
by 



hm if (r, fa 



Airr 



d9 In 



l + exp(-e,(0,r, fa 1 )) 
l + exp(-Ej(0,r, //£)) 



ri[e(0)|?H0)p 
d9 



(3.23) 
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after a partial integration. For simplicity we also assumed here parity invariance, that is 
|T"(0)| = \T£*(6)\. The derivation of the analogue to ( 3.23| ) for the situation when parity is 



broken is of course similar. Taking now the potentials at the end of the wire to be p^ 
V/2, the conductance reads in this approximation 



limG«(r) 

r—>0 



Qi 

2irr 



de 



dsi(e,r,V/2) 



l + exp[ £i (0,r,O)] dV 



v=o 



d[e(e)\T?(e)f 
de 



(3.24) 



In order to evaluate these expressions further, we need to know explicitly the precise form of 
the transmission matrix, i.e. the concrete form of the defect. An interesting situation occurs 
when the defect is transparent or rapidity independent, that is |T"(#)| — > \Tf\, in which case 
we can pursue the analysis further. Noting that de(9)/d9 = 25(9), we obtain 

Hi \T?\ 2 



lim Gf 



(r) 



de l (0,r,V/2) 



v=o 



Trr 1 + exp£j(0, r, 0) dV 
The derivative (fei(0, r, V/2)/dV can be obtained by solving recursively 

<fei(0,r,V/2) r ^ AT 1 dej(0,r,V/2) 



dV 



l + expej(0,r,V/2)] dV 



(3.25) 



(3.26) 



which results form a computation similar to a standard one in this context Q leading to the 
so-called constant TBA-equations. Here only the asymptotic phases of the scattering matrix 
enter via iVjj = lim^oofln^ij^— e)/Sij(e)]]/2iri. The values of £j(0, r, 0) needed in ( |3.25| ) can 
be obtained for small r in the usual way from the standard constant TBA-equations. 



3.5 Free Fermion with defects 

Let us exemplify the general formulae once more with the free Fermion. First of all we note 
that in this case in the TBA-equations ( |3.18| ) the kernel tp,-^ is vanishing and the equation is 
simply solved by 

£i(9, r, pLj) = rrrii cosh 9 — rp i . (3.27) 
Therefore, we have explicit functions for the densities with ( 3.19| ) and ( |3.17j ) 



1 n r ,n \ m i COsh 9/2-K 

Pi{0,r, Pi) = —mi cosh 9 and p i (9,r,p i ) = — — -. (3.28) 

2-7T 1 + exp(rmj cosh 9 — rp t ) 

According to ([O]) the direct current reads 

oo 

I«(r,V) = | Jd9[p^9,r,V/2)\T? { 9)f-p:(9,r,-V/2)\T? { 9)f 

— oo 

-pl(e,r,-V/2)\f?(e)\ 2 +pl(e,r,V/2)\ff*(0)\ 2 ] . (3.29) 
Using atomic units m e = e = fi = mi = qi = 1, we obtain explicitly with fl3.28p 

oo 

I«(r V) — — [ de ™™^HrV/2)\T<*(9)\* 

v > ' / ttJ P cosh(rcosh0) + cosh(ry/2)' { ' 

o 



20 



for |Tf (0) 
results to 



\T? {9) | = \ f{* (9) | = | (0) | = \T a (9) | . Then by Q the conductance 



G"(r) 



nu- 



ll 



d0 



coshfl |r tt (fl)| 2 

1 + cosh(rm cosh 9) 



(3.31) 



in this case. We have re-introduced dimensional quantities instead of atomic units to be able 
to match with some standard results from the literature. The most characteristic features can 
actually be captured when we carry out the massless limit as indicated in section 2.3.2, which 
can be done even analytically. Substituting t = e® , we obtain 



lim G a (r) 

m— >0 



h 



dt 



\ r TOL 

\ 1 L/R 



(ty/r)f 



1 + cosh(t) 



-L/R 
pa 
L/R 



(ty/rW 



\T?,Jy/r 



for y ^> r 
for y <C r 



(3.32) 



We have identified here two distinct regions. When y <C r we can replace the left/right 
transmission amplitudes by their values at y/r = 0. When y>r the transmission amplitudes 
enter the expression as a strongly oscillatory function in which y/r plays the role of the 
frequency. It is then a good approximation to replace this function by its means value as 
indicated by the overbar. It is straightforward to extend the expression 

( gjp to the 

case 

when the assumption on T a in ( |3.30| ) is relaxed and to the case with different values of y. To 
proceed further we need to specify the defect. 



3.5.1 Energy insensitive defects, T> (ijj,ip) = 0, T>P(i{>,ip), V(ip,iJ}), T> 5 ±(Tp,ip) 
Let us first consider the easiest example, which supports the general working of the method. 
When the defect is transparent, i.e., \T a \ = 1, we can compute the expression for the con- 
ductance (3.31) directly in the large temperature limit and obtain the well known behaviour 



lim 

+0,|T a |- 



G°(r) 



h 



[I 



(3.33) 



Alternatively we obtain the expression ( |3.33 ) also from equation ( 3.25|) and (3.27). In the 
massless limit of ( |3.32| ) we obtain e 2 /h which coincides with the result in [13|. However, we 
should stress that we consider here purely massive cases and the massless limit only serves as 
a benchmark. Note that a transparent defect in this context does not necessarily mean the 
absence of the defect. Considering for instance the defect P^(^/>,V ; ), we compute with ( |2.48| ) 



and (2.49) the same conductance as if there was no defect at all. Similarly simple are the 
computations for the defects V 1 (^ , ip) , T> s± tp). We simply get 



G°(r) = G^{r) = G 7 (r)cos 2 £ = G s± (r)/(1 + 4tan 2 B/2) 



h 



(3.34) 



Since the amplitudes do not depend on the rapidities, the TBA-kernel is zero and there is no 
contribution from this defect to the free energy, even unit length. 



3.5.2 The energy operator defect V a (tjj , tp) = gipifj 

For this defect the computation of the conductance according to ( |3.31| ) is more involved. The 
results of our numerical analysis of the expression (|3.31 ) are depicted in figure 4. 
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Figure 4: Conductance G(r) for the complex free Fermion with the energy operator defects as a 
function of the inverse temperature r, for fixed effective coupling constant B and (a) for varying 
amounts of defects I = 0,1, 2, 4. (b) for £ = 2 for varying distances y. 

We observe several distinct features. First of all it is naturally to be expected that when 
we increase the number of defects the resistance will grow. This is confirmed, as for fixed 
temperature and increasing number of defects, the conductance decreases. Second we see 
several well extended plateaux. They can be reproduced with the analytical expressions 
obtained in the massless limit ( 3.32j) . To be able to compare with ( 3.31 ) we re- introduce 



atomic units for convenience, i.e. e 2 /h — > 1/2-k. For a single defect there is only one plateaux 



and from (3.32) we obtain with fl2.43| ) 



G a (r) 



cos 2 B 
2vr 



(3.35) 



For B = 0.5 the value 0.1226 is well reproduced in figure 4(a). The lower lying plateaux 
correspond to the region when y <C r. In that case we obtain from ( 3.32| ) together with the 
expressions ( |2.44 )-(2.47) for a double and four defects 

2 

for y <C r , (3.36) 



G aia2 (r) 



1 

2^ 
1 

2^ 



cos 2 B 



1 + s'm z B 



cos 4 B 



cos 4 5 



for y <C r. 



(3.37) 
2 and 



2(1 + sin 2 B) 2 

For B = 0.5 the values 0.0624 and 0.0095 are well reproduced in figure 4(a) for 
£ = 4, respectively. The plateaux extending to the ultraviolet regime result from ( 3.32| ) and 
by taking in (|2.44) - (|2.47|) the mean values 

1 + sin 4 B 



G aia2 {r) 



7T (COS 2 (25) 
1 

+ 



3) 2 ' 

cos 8 B 



4vr 4tt[cos 4 5-2(1 + sin 2 B) 2 ] 2 ' 



for y S> r , 
for y 3> r. 



(3.38) 
(3.39) 



Also in this case the values for B = 0.5, i.e., 0.110784 and 0.084311 for £ = 2 and i = 4, 
respectively, match very well with the numerical analysis. Finally we have to explain the 
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reason for the increase from one to the next plateaux and why the curves are shifted precisely 
in the way as indicated in figure 4(b) when we change the distance between the defects. 
This phenomenon is attributed to resonances as we shall discuss in more detail in the next 
subsection. 



3.5.3 Resonances versus unstable particles 

In El we demonstrated that resonances may be described similarly as unstable particles. 
The latter provide an intuitively very clear picture which explains the relatively sharp onset 
of the conductance with increasing temperature. The temperature at which this onset occurs, 
say Tq can be attributed directly to the energy scale at which the unstable particle is formed, 



since then it starts to participate in the conducting process. The Breit-Wigner formula [47] 
constitutes in this case a relation for the mass M and the decay width T of an unstable particle. 
Supposing that in the scattering process between particles of type i and j an unstable particle 
can be formed, this is reflected by a pole in Sij(6) at = a — id. Then, for large values of the 
resonance parameter a one may approximate 

M 2 ~ l/2mjm,(l + cos<r) exp \a\ and T 2 ~ 2mjmj(l — coscr) exp |<r| . (3.40) 

Since a renormalization group flow is provided by M — > rM, one should observe that the 
quantities M ~ rie ai ^ 2 = r2e CT2//2 and T ~ rie CT1//2 = r2e CT2//2 remain invariant. Accordingly, 
this creation of the unstable particle should be reflected in the conductance as 

G(ri,a 1 ) = G(r 2 ,tr 2 ) for ne ffl/2 = r/ 2/2 . (3.41) 

This means we can control the position of the onset in the conductance by M and its extension 
in the temperature direction by T. For a model which possesses unstable particles we found 
indeed such a behaviour |46[| . From the data of the previous subsections we find that the 
conductance scales as G(rx,y{) = G{r2,y2) for r2yi = r\yi- Then the comparison with ( |3.4l|) 
suggests that we can relate the distance between the two defects to the resonance parameter 
as a = 2 ln(const/y). From the maxima in |T(0)| we may identify various as and in fact in 
this case the net result can be attributed to two resonances [E3] . 



3.5.4 Multiple plateaux 

Up to now, we have observed that we always obtain essentially two plateaux in the conduc- 
tance, no matter how many (> 2) and what type of defects we implement. The natural 
question arising at this point is whether it is possible to have a set up which leads to a more 
involved plateaux structure? It is clear that if we had many defects in a row separated far 
enough from each other such that the relaxation time of the passing particles is so large that 
they could be treated as single rather than multiple defects, then any desired type of multiple 
plateau structure could be obtained. In this case the conductance is simply the sum of the ex- 
pressions one has for each defect independently. Recalling the origin of the different plateaux, 
there is another slightly less obvious option. The density distribution function p r is a peaked 
function of the rapidity and if the resonances in T a (6>) would be separated far enough, such 
that they are resolved by p r , we would also get a multiple plateaux pattern. However, tuning 
the distance between the defects or the coupling constant will merely translate the position 
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of the resonances in the rapidity variable or change their amplitudes, respectively (see section 
2). 
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Figure 5: Conductance G(r2) for the complex free Fermion with the energy operator defects as a 
function of the inverse temperature r2, for fixed effective coupling constant B = 0.5 and varying 
temperature ratios in the two halves of the wire. 

Therefore the last option left is to change the p r s, which is possible by varying the temperature. 
Choosing now a configuration as in figure 3 with different temperatures T\ and T2, one can 
"create" a second plateau at half the height of the original one. The reason for this is simply 
that the cooled half of the wire will cease to contribute to the conductance as can be directly 
deduced from ( p. 31 ). We depict the results of our computations in figure 5. 
From this it also obvious that if we only cool the fraction x of the wire, the lowest plateau 
will be positioned at the height x times the height of the upper plateau. Thus, by combining 
these different configurations, i.e., different temperatures or defects, we could produce any 
desired plateau structure. 



4 Conductance from the Kubo formula 

Having computed the DC conductance by means of a TBA analysis, we want to proceed now 
by introducing an alternative method for the acquisition of the same quantity, that is the 
evaluation of the celebrated Kubo formula]]] Q 

00 

G(T) = - lim ^ / dt ^ ( J ( f ) J (°))T, m ■ (4-1) 



The key quantity needed for the explicit computation of (4.1) is the occurrence of the current- 
current correlation function (J(r) J(0)) Tm . In the latter, the subscripts (T, m) indicate that, 
in general, one is interested in a situation when both, the mass scale of the particles in the 
quantum field theory and the temperature, are non-vanishing. This is precisely the same 
regime in which we have carried out the TBA analysis in the previous section and ultimately 



"'"For a model independent derivation in the context of dynamical respose theory see, e.g., 
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we want to compare the outcome of both computations. So far, formula (4.1) still refers to 
a situation in which no defect is present in the theory. Later on we will see how the Kubo 
formula can also be generalized in order to incorporate the presence of defects. 

As a consequence of the central role played by the two-point function of the current 
operator in (4.1), we will devote an important part of this section to recall the key features 



of a concrete method which will allow for the computation of such a quantity that is, the 
form factor bootstrap approach |Tlj| . To carry out this program one needs essentially as 
the only input the scattering matrix and it is then in principle possible to compute form 
factors associated to various local operators of the quantum field theory under investigation. 
Form factors are defined as matrix elements of some local operator 0{x) located at the origin 
between a multiparticle in-state and the vacuum, 

fS W -^{6 1 ,9 2 . . . , n ) := (010(0)14^)4^) . . . 4 n (0 n )>- (4-2) 
They can be obtained by a direct computation once a representation for the operator involved 



is known or as solution to a certain set of physically motivated consistency equations [11, H 



48, HJJ, in a similar fashion as one can determine exact scattering matrices or transmission 



and reflection amplitudes for 1+1 dimensional integrable systems as discussed in section 2. 



In the zero-temperature regime, the latter fact is well-known since the original works [11] 
and has lead successfully to the computation of correlation functions for many models, albeit 
in most cases only approximately. It is easy to show that once the corresponding form factors 
associated to two local operators O and O' are known, the computation of their two-point 
function is reduced to the task, still non-trivial, of evaluating the following series 



'1 " " " du n -|-|- ^-rrrii cosh 0, 



<o(>W0)> Trf . m = £ £ /%^IL 

xF OW..a n{e ^ A) ^O'W..^^ * (43) 

with = (— ir, 0). The previous expression is simply obtained by introducing a sum over 
a complete set of states in between the two operators involved in the correlation function 
and by shifting the operator 0(r) to the origin thereafter. However, as indicated by the 
subscripts, this formula applies only to the zero temperature regime. Obviously, when setting 
O = O' = J, the correlation function Q4.3|) is precisely the quantity entering the Kubo formula 



(4.1), although for T = 0. It is therefore necessary to find a generalization of the expansion 
( |4.3| ) to the (T / 0)-regime. Such type of generalization was first suggested in |5(]] for the 
Ising model. It appears, however, to be difficult to generalize this to dynamically interacting 



models [51| and since by now this has not been achieved we shall concentrate on the zero 



temperature regime in this paper. 

4.1 Conductance through an impurity 

With the help of (^^) we could in principle compute the current-current correlation function 
and therefore evaluate the Kubo formula when there are no boundaries or defects present. 
With regard to the inclusion of boundaries, the first examples in which the Kubo formula 
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was generalized in order to accommodate that situation were provided in fig] . In there the 
expression (|4.1| ) was evaluated for the sinh-Gordon and sine-Gordon model in the m = T = 
limit and in the presence of a boundary. This was done by replacing the vacuum state |0) 
with a boundary state \B) in the current-current correlation function as follows 

(J(r) J(P)) T ,m-+ ( J ( r ) J (°) B )T,m = <0| Ar) J(0)B |0) Tim . (4.4) 

The boundary state \B) :=B|0) is understood as the action of a boundary operator B on the 
vacuum state |0) ^9|. Following p^j, one exchanges usually the roles of space and time, such 
that the correlation functions are radially rather than time ordered. This is the reason why 
the boundary operator B can only enter at the very right or left, since one formulates such 
theories in half space. In contrast to the boundary a defect can also enter in-between the 
operators. Therefore, in order to include a defect in ([4.1|), one has to consider terms of the 
form 

(J(r)J(0)) T:m ^(J(r)Z a J(0)) Tm , (J(r)J(0)Z a ) Tm , (Z a J(r) J(0)) T>m , (4.5) 

where Z a represents the defect operator. 

As a consequence of (4.5), the evaluation of the defect Kubo formula will require the 



computation of matrix elements involving the operators Z a 

oo 

G a (T) = -lim^ / dr e^ r (J(r) Z ai ■ ■ ■ Z an J(0)) T , m . (4.6) 

— oo 

Equation ( ]Q| ) expresses the conductance for a situation in which n generic defects Z ai ■ ■ ■ Z an 
are present in the theory and located at positions y ai ■ ■ ■ y an in space. The defect degrees of 
freedom are encoded into the vector a ={«i, • • • , a n }, as done in previous sections. In order 
to compare with the TBA results we would like, of course, to compute the conductance in the 
massive, finite temperature regime. As mentioned the evaluation of temperature dependent 
correlation functions is still poorly understood, even for the simplest models. In addition, 
the presence of the limit in the parameter oj, together with the introduction of the defect 
operator Z a in the current-current two-point function makes the generic evaluation of ( |4.6|) 
fairly involved and constitutes a problem which in general can not be solved analytically. This 
is specially cumbersome when double defects are considered, since the expressions for the 
reflection and transmission amplitudes ( 2. 21] ), (|2.22 ) are, in general, quite messy to handle. 



For these reasons it is interesting to start with a more simplified situation, in which some 
analytical calculations can still be performed, that is the T = regime. One may now view 
( ]4.5| ) as a three-point function and extend the expansion fl4.3|) to the case when three operators 
enter the correlation function. This will only require the inclusion of one more set of complete 
states, such that (]4.5|) is expanded in terms of the form factors of the three operators involved 



n,m=l Pi—Pn 
v\—v m 

x (z^n) . . . z H {e x )\z a \z vl Cei) . . . z Vm (e m )) f^- v ™(o x ,. . .,~e m )* ( 



-r m i cosh 6i 
■ e i=l 



26 



We will now restrict ourselves further and consider the massless version of (4.7). In this limit 



the results obtained for the conductance should agree with the UV-limit of the conductance 
computed by means of (|3.1|), Q3.31 ). Such a limit can be carried out by exploiting the mass- 



less prescription suggested originally in |37[ and already introduced in the paragraph before 
equation ( |2.43| ). For the form factors in (|4.7| ) the massless limit yields 

lim F^---^(e 1 + Vl a,...,e n + Vn a) = F^ 1 h-f-(e 1 ,...,e n ), (4.8) 

a— >oo 

with rj i = ±1 and V{ = R for r] i = + and V{ = L for rj i = — . Namely, in the massless limit 
every massive n-particle form factor is mapped into 2 n massless form factors. Using these 
expressions, performing a Wick rotation and introducing the variable E = X)ILi > we 



obtain from (4/7) 



n,m=l Ml— Mn ' ' V ' 

V\—V m 

x (z* (e B ) . . . z* x {e x )\z a \z*(e x ) . . . z*je m )) f^j^, . . . , e n )*e+ B . 

We note that for the massless prescription to work, the matrix element involving the defect 
Z a can only depend on the rapidity differences, which will indeed be the case as we see below. 
Performing the variable transformation 6 n — * In E'/rh n — ^™ = iThi/rh n e 9i , we re- write the 
r.h.s. of @ as 

^ In E'/rhn oo _ 

E E I * J / ^^S^^-Mmm 

x (zfoe n {Ef)) . . . ft) . . . z?je m )) J^?j£»(0i, • • • , U*^'. 

We substitute now this correlation function into the Kubo formula, shift all rapidities as 
G% Oi+ ln£"/m n , Q% — * 9i+ ln£"/m n use the Lorentz invariance of the form factors^] and 
carry out the integration in dE' 

oo _ 

w 2 S -2 r d0f- dQ n -i f d0 a • • ■ d9 m 1 



G « = _ lim " y 



^ J n!(2vr)« 7 m!(2vr)- l-JXl^i/^* 

,7 1 I ...,7 n — oo — oo 

Z*(ln(l - Y^'l ^M"^)) • • • Z*(Pi)\Z a \Z*(jh) ■ ■ ■ Z? m {O m )) (4.11) 



n-l 



^Denoting by s the Lorentz spin of the operator O and A being a constant, the form factors satisfy 
F e>K...M„ (ei + A, . . . , # n + A) = e sA F^^-^^i, . . . , n ) . 
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We state various observations: Since the matrix element involving the defect only depends 
on the rapidity difference, it is not affected by the shifts. The Lorentz spin s = 1 plays a 
very special role i n ((4.11| ), which makes the current operator especially distinguished. In that 
case the r.h.s. of ( }4.1l| ) becomes independent of the frequency u and the limit is carried out 
trivially. Furthermore, since the final expression has to be independent of rh n , we deduce that 
the form factors have to be linearly dependent on rh n . 

One may now compute the form factors by solving either the associated consistency equa- 
tions or by using concrete realizations of the operators. For those form factors involving the 



current operator J, a realization in terms of the ZF-algebra was given in [28] for complex 
free Fermion type models and used to compute the corresponding matrix elements. We will 
determine form factors involving the defect operator in the same fashion, which means we 
require a concrete realization for the operator Z a . 

4.2 Realization of the defect operator 

A realization of Z a can be achieved very much in analogy to a realization of local operators, 



i.e. as exponentials of bilinears in the ZF-operators [52|. For the case of a boundary a generic 



model independent realization for the boundary operator B was originally proposed in [29] for 



the parity invariant case, i.e., R = R . This proposal was generalized to the defect operator 



in [32] with the same restriction and for self-conjugated particles. This realization was used 
by the authors for the computation of various matrix elements involving the defect operator. 
Our aim in this section is to extend this realization in order to incorporate the possibility of 
parity breaking as well as non self-conjugated particles. A non-trivial consistency check for 
the validity of our proposal will be ultimately provided when exploiting it in the computation 
of the conductance, obtained before by entirely different means, that is the TBA approach. 



The realization we want to propose here is a direct generalization of the one presented in [32], 
namely 

1 f 00 

Z a =:exp[— / D a (0)d9] : , (4.12) 
where : : denotes normal ordering and the operator D a {0) has the form 

D a {9) = Y^[K?(e)z}(e)zU-9) + Km*Zi(-m(.o) 

i 

+w l a (e)z}(e)z i (9) + wneyzk-ew-o)] , (4.13) 

with 

7 77- _ _ 777- 

K?(0) := Rf{- - 9), K?(9) := Rf(- - 9), (4.14) 

W?(9) := T?(- - 9), Wm := T?{- - 9). (4.15) 

In comparison with [^] we have used a slightly different normalization factor, since in general 
we have contributions in the sum over i in ( 4.13j ) including both particles and anti-particles, as 



for the complex free Fermion we shall treat below. Following the arguments given in [29|, the 



operator D a (9) depends on the amplitudes R(0), T(6), R{9) and T{9) with their arguments 
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shifted according to ( (4.14D - fl4.15|) , as considered also in [^5], The reason for these shifts is 
the exchange of the roles played by the space and time coordinates = (t, x) — ► (ix, it), which 
was already mentioned after (4.4). Doing this and keeping simultaneously the product x^ 



invariant requires the rapidity shifts in (|4.14|) -( |4.15 ). In our context, this implies that we must 
now not only perform the shifts (4.14)-(4.15) in our expressions, but also, with regard to the 
positions of the defects, the change y a — > iy a should be implemented. The latter replacement 
will play an important role since, similar as in the TBA case, the amplitudes (p^)-( |4.15Q 
will become in this way strongly oscillating functions of y a . Therefore, we may be able to 
carry out once more certain analytical calculations, by replacing the mentioned functions with 
their mean values. In ( 4.13; ) we have already specialized to the case when the reflection and 
transmission amplitudes are diagonal both with respect to the particle and defect degrees of 
freedom, since that will be the situation for all the examples we want to treat in this paper. 



4.3 Defect matrix elements 

Having now a concrete generic realization of the defect ( 4,13j ), we can compute the defect 
matrix elements. One way of doing this is to solve a set of consistency equations which 
relate the lower particle matrix elements to higher particle ones, similar as in the standard 



form factor program plj| . Such kind of iterative equations were proposed in [25] for a parity 
invariant defect and for a real free fermionic and bosonic theory. First we note that the 



operator ( 4.12 ) becomes 



in the limit R = R = and T = T 
operator and, according to fl2.2|) , 



1 f°° 

z a =-. ex P [— J de^2 z}(e)z i (e) ] ■., (4.ie) 

1. The defect should act in this the identity 



(Z i (0 1 )Z a Zl(9 2 )) = 2ir5(6 12 )5 



ij- 



(4.17) 



holds, simply by employing Wick's theorem when carrying out the necessary contractions. 
For two particles we find, 



(^(0i)Zi(0 2 )Z a ) 

(z a zUe 1 )zHe2)) 



(Z i (9 1 )Z a z](9 2 )) 



ttK?(9 2 )6(9 12 ), 

7tK^(9 1 )*5(9 12 ), 

nW^{9 l )5{9 12 )5 iv 



where we recall from section 2 the notation 6i 2 = 
convenience we have introduced the functions 

K«(0) = K^{e) + S u (-29)KS*(-9) 
W?(6) = W<*{9) + W?{-9)* = W?{ 



01 + 9 2 and 9 
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(4.18) 
(4.19) 
(4.20) 

?i — 9 2 . For later 



-K?(e) + S*(2e)K?(-9), 
-9) +W T a (9)* = W-?(9)\ 



(4.21) 
(4.22) 



since the K", K", and Wf amplitudes defined by ( |4.14 )-(4.15) will repeatedly appear in 
the combinations ( 4.21| ), ( |4.22| ) in what follows. The latter equalities in ( 4.21 ), ( 4.22| ) follow 
simply from 



W?(9) = W-?{-9) = Wf(in - 9)*, K?(9) = S a (29)K?(-9) = S«(29)K?(iir 



(4.23) 
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which are in turn consequences of the crossing-hermiticity properties (2.16)-( p.l7|) . Having 
these matrix elements we can construct the ones involving more particles recursively from 



F, 



H m ...fJ. 1 U 1 ..M n 



9t,9[--- O := (Z„j9 m ) . . . Z^{9 x )Z a . . . Zt n (9' n ) 



1-1 



w(m<(0i) n wm^< 



ct 



(9 m ...9 l ...9 2 ,9[...9' n ) (4.24) 



1=2 



p=l 



l-l 



+^E<W(#i - ^X(^i) II S^ p {9 lp )F^-^- 

p=l 



(9 m ... i 



1=1 



'2,Cl 



...9i,9[... 9 n ) 



(4.25) 



i-i 



(# m ... I 



'1,^2 



Z=2 



p=l 



Z-l 



■^J^M " 9i)W^(9[r J] S VWp (0 lp )F£ 
i=i p=i 



m ...H l ...(l 1 i/2...Vn, 



h,9 2 . ..9' n 



Here we denoted with the check on the rapidities 9 the absence of the corresponding particle 
in the matrix element. It is clear from the expressions ( [4.12 ) and ( 4.13 ) that the only possible 
non-vanishing matrix elements ( [4.24 ) are those when n + m is even. Taking (4.18)-( 4.20| ) as 
the initial conditions for the recursive equation ( 4.24[ )-( 4.25 ), we can now either solve them 
iteratively or use (4.12) and evaluate the matrix elements directly. 



4.4 Free Fermion with defects 

Similar as for the TBA we want to exemplify our general formulae with the complex free 
Fermion. We consider now the particularization of the defect realization fl4.13| ) to this case. 
Then the generators of the ZF-algebra Zi(9), z\{9) are just the usual creation and annihilation 
operators Oi(9), a\(9) in the free fermionic Fock space and we have to distinguish between 
particles i and antiparticles i. For the complex free Fermion it is interesting to notice that 
the realization ( 4.12 ) resembles very much the one employed in [52, ^] for a prototype local 
field. 



4.4.1 Defect matrix elements 

Let us now use fl4.12| )-( [4.13| ) in order to evaluate matrix elements involving the defect operator. 
In what follows, the most relevant matrix elements are those involving four particles, for which 
we compute 

(a l (9 1 )a- l (9 2 )Z a at(9 3 )al(9 4 )} = w^ l {9 l 9 2 )5{9 lA )5{9 23 ) + A£(0i,0 4 )<KM<*(M, (4.26) 
(a l (9 1 )a l (9 2 )Z a a](9 3 )a](9 4 )) = -^ 2 W^{9 l )W^{9 2 )8{9 l3 )8{9 24 )5 ij , (4.27) 



{a l {9 1 )a k {9 2 )a i {9 3 )Z a a\{9 i )) = K 2 W?(9±)k?{-9 2 ) [5(9 U )5(9 23 ) - 5(9 12 )8(9 U ) 
(ai(9i)Z a al(9 2 )al(9 3 )4(94)) = n 2 W^{9x)k^{-9 z y \5(9 23 )5(9 14 ) - 5(9 12 )5(9 U ) 
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with the abbreviations 

<(0i,02) = 7T 2 wr(e 1 )wr(e 2 ) ffiVM = i?k?{e x )k?{o % y . (4.28) 

One can also find solutions for all n-particle form factors either from Q4.24[ )-( fi.25p or by direct 
computation. For instance we compute 

. . min(n,m) ^ m +n-2k ^n+m roo 
ptnx jj nx ti ( n n n' a> \ — \ ' / AR AR 

F a (9 2m ...9 1 ,9 1 ...9 2n )- ^ - _ _ _ _ J_J^...d0 2n+2m 

x det^l 2n (/3 1 . . . P 2n ; 9[ . . . 9' 2n ) detA 2m (P 2n+1 . . . P 2n+2m \ 9\ ■ ■ ■ 9 2m ) 
k 

x II Wri^WriP^-M^ - P 2n+2p W 2p _ x - 02n+2 P -i) (4-29) 
p=i 

n n+m 

x II K^(P 2p )*6(p 2p + p 2p _ 1 ) I] K?(p 2p )5{f3 2p + I3 2p ^), 

p=l+fc p=l+fe+n 

where -4^(#i . . . Of, 9 X . . . 9'^) is a rank £ matrix whose entries are given by 

A% = co S 2 [(i - j)7r/2]5(9 i l<i,j<£. (4.30) 



The matrix elements are computed similarly as in [£8| and references therein. Likewise we 
compute 

-n( i \n— 1 roo n 



n(_-\\n-l roo " ^ 

jF nx l+m x i( ^ . . . g ^ ) = ,5 n , m ^-^ / dft . . . d(5 n J] 

n - J-oo , , 



/,-.! 



x det £ n (0 n . . . i; 1 . . . fi n ) det B n {fi x . . . /3 n ; 9' x . . . 9' n ), (4.31) 

where we introduced a new rank i matrix B i (9 x . . . 9g\ 9 X . . . 9'g) whose entries are now simply 
given by 

B £ ij = S(9 i -9' j ), l<i,j<£. (4.32) 

Since ( |4.31| ) is simpler than ( (4.29|) we use it to demonstrate explicitly that it satisfies the 
recurrence relations (|4.24| ) and ( [4.251) . The other cases work the same way. Starting with the 
expansion of the determinant det,6 ra (# n . . . 9 X ; P x ■ ■ ■ (3 n ) with respect to the row involving the 
variable 9 X gives 



det B n (9 n ...0 X ;(3 X ... P n ) = Y J (-^) n+l+l K0i-Pi) det B n ~\9 n . . . 9 2 ; (3 X . . . ft . . . n ). (4.33) 

i=i 

Inserting then ( |4.33| ) into ( f4.31| ), we obtain 

n n roo 

F nxi + mxi {9n du ^ e ,j = ^ £ wr{di) / ^ ^ 0n 

n. J=1 J-oo 

l-l n 

x(-l)'IIW ]I W?(O k )detB n - 1 (O n ...9 2 ;P 1 ...0 l ...p n ) 

k=l k=l+l 

x det B n (f3 x ...p l ^9 l ...p n ;9' l ...9' n ). (4.34) 
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Expanding now the second determinant in (4.34) with respect to the l-th row, which involves 



the rapidity 9\, and using the fact that the (3s are just integration variables and therefore, the 
sum in I gives actually n times the same contribution, we can write 

n n roo 
F nxi+mxi {9n ...0 X J 1 ...# m ) = 5n, mr —^ £ WfW^l ~ 9' p ) / df3 l ... d{3 n _ x 

\ n )■ x J-oo 



P= 

p— 1 n 



x(-l) p n^°W II W l a (6 k )detB n - 1 (6 n ...e 2 ;P 1 ...(3 n _ l ) 

k=l k=p+l 

x det B n -\^ . . . $ n _xA ...e' p ...0' n ). (4.35) 

We recognize now the matrix element with two particles less on the l.h.s. of ( |4,35| ) and can 
re-write it as 

n 

F™ i+n * i (d n ...O l ,e' 1 ---0'n) = ^(-If^CWfli-flp) (4.36) 

p=l 

Xfa (?n-l ■ ■ ■ V2,Vl ■ ■ ■ Vp ■ ■ ■ V n )i 



which is in complete agreement with ( [4.24 ). The validity of (|4,25 ) can be checked similarly. 



4.4.2 Conductance in the T = m = regime 

It is well-known that for a free Fermion theory (also for a single complex free Fermion) the 
conformal [/(l)-current-current correlation function is simply 

<J(r) J(0)) T=m=0 = 1. (4.37) 

This expression can also be obtained by using the expansion (|4.3j), together with the massless 
prescription outlined before ( 2.43j ) (see [pi] ) and the expressions for the only non- vanishing 



form factors of the current operator in the complex free Fermion theory 

F^\6, 0) = -F 2 J{i \6, 0) = -inme^ . (4.38) 

In particular, the massless limit of the previous expressions gives, according to the massless 
prescription, 

F$g(e,0) = -F$g(9,0) = -2*ime t P, (4.39) 

F$(0,0) = F^(9,9)=F^\9,9) = 0, (4.40) 

F J L f{e,~e) = F^(9,9)=F^(9,9) = 0. (4.41) 



Inserting ( [4.37 ) into ( |4.1[) reduces the problem of finding the Fourier transform of the function 



2 which is given by V dr e tujr r 2 = —iru for uj > 0, with V denoting the principle value. 
This yields in the absence of a defect G°(0) = l/2n, in complete agreement with the limit 
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Let us now consider a more complicated situation, that is, the evaluation of (4.6) for 
T = m = in the presence of n defects Z ai ■ ■ ■ Z an located at positions y ai . . . y an in space. 
The correlation function ( |4.5[ ) can now be obtained with the help of (^^), which has to be 
generalized for three-point functions. This requires the inclusion of one more sum over a 
complete set of states in ( |4.3| ). Fortunately the only non- vanishing form factors of the current 
are (4.39), which means the expansion ( [4.3[ ) will already terminate for two particles. Explicitly, 
we find 



(J(r)Z ai ---Z an J(0)) 



T=m=0 



E 



d6id6 2 d0 3 d84, „j\ii to a , 

2(27r) 4 FrrViM 



b RR (#3,04, 



xe -rHe 6 ^2) {ai{0i)a _ MZai . . . z an 4(e 3 )4(e A )) m=0 , 



(4.42) 



In the light of the expressions ( 4.28| ), we can re- write (4.42) in a more explicit form without 
the need of specifying a concrete defect yet. Inserting ( 4.2C ) and fl4.39|) into ( 4.42 ), we find 



{J(r)Z ai ...Z an J(0)) 



T=m=0 



+ 



dOi 



m 
~2~ 



E 



d,e 



\ c —2rm cosh 6»i j£<x\R 



d,e 



2 j^ol\R 



(4.43) 



where we have exploited the crossing relations stated in (4.23). Here the functions 
K^ R {6), . . . defined in (4.43) are the massless limits of the corresponding functions W^{6), 
Kf c {6), . . . For all the defects we will consider below, it turns out that the first contribution to 
the previous correlation function is actually vanishing, so that ( [4.43| ) is considerably simplified. 
In many of the examples we will treat later, this is due to the fact that the amplitudes 
Kf-{9) are vanishing in the first place, as a consequence of the crossing relations ( 4.23[) . This 
will be the case for all energy insensitive defects for which we will present a case-by-case 
computation of the conductance below. The vanishing of the reflection part in ( 4.43|) also 
occurs in some cases as a consequence of the parity of the function Kf L {6). For instance, 
we find that, for the energy operator defect such function, although initially non- vanishing, 
satisfies k?{0) = -K?{-9), such that lim m _ S-oo d9 K^(d)* = 0. 

We can now either use Q4.43D in (4.6) to compute the conductance or evaluate the expression 
(4.11) directly in which the frequency limit is already taken 



G a (0) 



2(2vr) 



3 



d6e e w$ RR [ln(l 



(4.44) 



There are, in addition, further generic results which can be obtained independently of the 
specific defect. We present them at this stage and will confirm their validity below by some 
specific examples. Specializing to the case in which all I defects are of the same type and 
equidistantly separated, i.e., y = y ai = ••• = y an . As in the TBA context ( [3.32 ), we can 
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identify two distinct regions 



w 



a\RR, 



7T 



a\R, 



W 



for finite y 
for y — > 



(4.45) 



where we used in addition ( 4.22Q . Supported by our explicit examples below, we find that for 



y — > in ( 4.45| ) the amplitudes W 
we have already argued above 



at\R 



k 



(0) become independent functions of the rapidity. As 
"'^(01,02) = °- ( 4 - 46 ) 



The two regions specified in ( [4.45 ) are in complete agreement with the regions identified in 
equation (|3.32 ), since we also consider here the massless limit. When exploiting ( |4.45[) , our 
explicit examples below yield the values of the conductance as those computed from (|3.32j ). In 
the regime y — > this is very apparent, since when IW"™'^ = const it becomes equal to 2|T"| 
and the conductance reduces just to the constant factor given by (4.45) times the value \/2ir 
obtained when defects are not present in the theory. The vanishing of the first contribution 
in ((4.431) is also quite suggestive with regard to the TBA results, since the conductance 
obtained in terms of thermodynamic quantities only involves the moduli of the transmission 
or the reflection amplitudes, but not both simultaneously and, in the light of the previous 
discussion, that seems to extend also to the form factor computation. The coincidence in the 



regime for finite value of y between the Kubo formula based on ( 4.45 ) and the results from 
the Landauer formula are less obvious and we support this by some explicit computations for 
several specific defects, similar as in subsection 3.5. 

4.4.3 Energy insensitive defects, V°(ip,ip) = 0, P /3 (V> , if)) , £> 7 (^,V0, ^WsVO 
A simple example to start with, which at the same time provides a first test of the working of 
the Kubo formula in the UV-limit is a transparent defect, i.e., purely transmitting. As shown 
in subsection 2.3.3, examples for this are the absence of a defect P ^,?/') = as well as the 
defect T>^(Tp,ip) = g'tp^tp, for which the associated reflection and transmission amplitudes 
are given by Q2.48|) and (|2.49| ). In this case the observation K"(9) = is of course trivial 
and therefore only the second integral in ( |4.43| ) is relevant. The situation in which there is 
no defect was already commented on in the paragraph after equation ( 4.41 ). We found in 
that case that the Kubo formula leads to entirely consistent results with regard to our TBA 
analysis, that is G°(0) = 1/27T in the massless limit. Considering now a theory with n defects 
of the type P^ i (^>,?/>), we find 



TV", 



(0) = 2e 



-inB 



K 



(0)=O, 



(4.47) 



simply by exploiting the expressions ( 2.4S| ) for a single defect and the formulae fl2.23|) and 



(2.24). From (4.47) it follows that 



W s (01)02) 



4vr 2 



(0 



1,V2) 



0, 



(4.48) 



just, as in the case in which the defect is absent. Therefore, we recover once more the value 

G ft-A( ) = -L (4.49) 

Z7T 
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for the conductance at zero mass and temperature. 

The next complication arises for defects whose reflection and transmission amplitudes are 
simultaneously non- vanishing, but at the same time are independent functions of the rapidity. 
As we have seen in section 3.5, those defects can be very easily handled in the context of a 
TBA calculation for the conductance, since the modulus of the transmission amplitudes is a 
constant, depending only on the defect coupling g. Therefore, the conductance is simply given 
by the constant \T a \ 2 times the value fl4.49| ). The vanishing of the function K^(9) can be 
established in those cases by exploiting the crossing properties listened above. Namely, from 
( 2.16 )-( p.l7|) we find Ki = K* = —K^, whenever the reflection amplitudes are independent 



of the rapidity, and therefore Kf- = Ki + K* = 2Re(-?Q). The latter quantity is zero for 
the defects D 7 and D s± treated in section 2.3.4, when setting y = 0, since the reflection 
amplitudes are purely complex quantities. Therefore, the Kubo formula computation leads to 
the same results as found in subsection 3.6.1, since we find 

w]-! RR (9 1 ,d 2 ) = 4vr 2 (l + 4tan 2 J B/2), k]j R R {9 X , 9 2 ) = 0, (4.50) 
4vr 2 



w ^ RR (e 1 ,e 2 ) = k^ RR (e 1 ,9 2 ) = o, (4.5i) 

COS 1 ) 



which yields 



2tt 2ir cos z B 



(4.52) 



As expected, for B = we recover once more the value ( fl~49P . 



4.4.4 The energy operator defect T>(ip, ift) = gipip 

Let us now treat the energy operator defect, the example which has been most extensively 
studied in our previous sections. Considering first a theory possessing a single defect of this 
type, we find 

Wm= 4 l ;f C0Sh ? R aad Km= 2iS[ " BSmh h l (4.53) 
1 v ' cosh 29 + cos 2B 1 w sin B - cosh 9 

Therefore, in this case the amplitude Kf L (9) is non- vanishing. However, we find that Kf L {9) = 
—Kf c {—9). This means that the integral of this function (or its complex conjugated) is 
vanishing. Consequently, only the transmission part contributes non- trivially t o (p~43| ). In 
order to evaluate ( 4.26] ) in the massless limit, we are interested in this limit of ( 4.53j ) which 



enters equation ( f4.43| ). We obtain 

w^ RR (9 1 9 2 ) = 47T 2 cos 2 B, (4.54) 



which, together with (4.39) leads to the result 



cos 2 £ _^ na(n , _ cos 2 B 
again in agreement with the corresponding result (|3,35|) from the Landauer formula. 



(J(r)Z a J(0)) T=m=0 = =^ G a (0) = (4.55) 
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Let us now proceed to the study of the conductance in the presence of a double defect. 
Again, we consider first the case T = m = and two defects of the energy operator type 
located at the origin and at a distance y from the origin, respectively. Expression Q4.42| ) 
again holds for that situation with n = 2. As explained in the paragraph after equation 
(2.22), the Greek indices in the defect operator encode also the space dependence. The 
reflection and transmission amplitudes are computed according to ( 2.21| ) and ( 2.22j ) with 
( 4.53 ). These functions can thereafter be substituted into equation ( [4.26D in order to determine 
the explicit form of the functions w ^} a2 ^ RR and ^ a2 \ RR [ n (4.26) for the double defect system, 
which depend now on the distance y between the defects and their expressions become very 
cumbersome. Once more, it is possible to show that the contribution to the conductance 
depending on £. aia2 l R ^ j s vanishing and therefore only the function w;^ 1Q2 '^ R will be of further 
interest to us. However, it is relatively easy to show that in the massless limit we find 



aia 2 |Ri?//j n \ 



47T 2 cos 4 B 



1 + cos(2rhye ei ) sin 2 B 
1 + 2 cos(2mye dl ) sin 2 B + sin 4 B 
1 + cos(2mye 02 ) sin 2 B) 
1 + 2 cos(2mye 8 ' 2 ) sin 2 B + sin 4 B 



(4.56) 



such that we obtain 



{J(r)Z ai Z a2 J(0)) T=m=0 
G aia2 (0) 



4 7r 2 ? ,2 

2 1 + sin 4 B 



4 [1 + sin 4 B] 
r 2 [cos 2 {2B) - 3] 



2 - 



tt[3-cos 2 (2 j B; 



l2 ' 



(4.57) 
(4.58) 



which precisely agrees with the corresponding result ( |3.38| ) obtained from the Landauer for- 
mula. The overbar denotes as before the mean value of the corresponding function. 

As explained above, we can also predict the precise position of the second plateau obtained 
within the TBA analysis given in equation ( |3.36| ). This is achieved by considering previously 
to the UV-limit, the limit when the distance between the defects y — > 0. By doing so we find 



which gives 



lim (J{r)Z ai Z a2 J(0)) T=m=0 



4tt 2 cos 4 B 
(1 + sin 2 Bf 



1 



cos 4 B 



lim G aia2 (0) 



r 2 (1 + sin 2 Bf 
1 cos 4 B 
2tt (1 + sin 2 B) 2 ' 



(4.59) 

(4.60) 
(4.61) 



in agreement with the value (3.36). 

Finally, in order to match all the results in subsection 3.6, we would like to address also 
the case £ = 4 in (|4.6[), that is, we consider now a complex free Fermion theory with four 
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by y. For the first region in ( 4.45 ), that is the UV- limit, we find 



equidistant defects of the type V a (ip,ip) = giptp. As usual, we denote their mutual distances 

(4.62) 



w 



/l(0l)/l(02 



{Wx)-mi)\W2)-W2)y 



with 



(5 - cos 2B) cos(4mye e ) + 2 cos(6mye e )) sin 2 B) - 128vr 2 cos 4 B(2 + (6 cos(2mye e 
1192 cos 2B - 348 cos 4B + 24 cos 6B - cos 8B - 995 - 256 sin 2 B cos(6mye 9 )) 
128 sin 2 #((17 - 12 cos 2B + cos AB) cos(2mye e ) - 4(cos 2B - 2) cos(4mye e ) . 



This expression appears somewhat messy, but when proceeding as indicated in fl4.45|) it will 
simplify considerably. Computing the mean value of this function we find 



(J(r) Za\ Z a2 ^ci'i ^0.4 

J(0)) T= 



m=0 



G 



RR 


1 


47r 2 r 




2^2 




1 + T 


JL( 


4?r \ 



(0i,e 2 



cos 8 £ 



[cos 4 B - 2(1 + sin 2 £) 2 ] 2 
cos 8 5 

4vr ' [cos 4 £ - 2(1 + sin 2 B) 2 } 2 



(4.63) 
(4.64) 



in complete agreement with the corresponding TBA value (|3.39| ). We can also predict the 
precise position of the second plateau which, according to (|4.45| ) is expected for the conduc- 
tance. Once more we find complete agreement with the outcome of our TBA analysis, since 
in this case 2 
hm< ia2Q3a4 (0!0 2 ) = ( 27rcos 4 £ \ 

y^o 11 v x ' z > Vcos 4 B - 2(1 + sin 2 B) 2 J ' v ; 



which gives 



lim (J {r)Z ai Z a2 Z a3 Z a4 J (0)) T=m=0 

y^O 



lim G aia2a3a4 (0) 

y^O 



1 / cos 4 5 

r 2 Vcos 4 £ - 2(1 + sin 2 B) 2 

1 / cos 4 5 

2^ Vcos 4 S - 2(1 + sin 2 B) 2 J ' 



(4.66) 
(4.67) 



that is, the same expression as (|3.3?1 ). 



5 Conclusions 

We have exploited the special features of 1+1 dimensional integrable quantum field theories 
in order to compute the DC conductance in an impurity system. For this purpose several non- 
perturbative techniques have been used. As the main tools we employed the thermodynamic 
Bethe ansatz in a Landauer transport theory computation and the form factor expansion in 
the Kubo formula. 
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The comparison between the Landauer formula and the Kubo formula (l.k) yields in 



•particular an identical plateau structure for the DC conductance in the ultraviolet limit. 

We have explained to what extend integrability can be exploited in order to determine the 
reflection and transmission amplitudes through a defect. Unfortunately, for the most interest- 
ing situation in this context, namely when R/R and T/T are simultaneously non- vanishing, 
the Yang-Baxter-bootstrap equations narrow down the possible bulk theories to those which 
possess rapidity independent scattering matrices ^6|. By means of a relativistic potential 
scattering theory we compute for several types of defects the R/Rs and T/Ts, thus enlarging 
the set of examples available at present. We confirm that for real potentials parity is pre- 
served, but otherwise essentially all possible combinations of parity breaking can occur. From 
the knowledge of the single defect amplitudes the multiple defect amplitudes, which exhibit 



the most interesting physical behaviours, can be computed in a standard fashion |33| , 34 1. 

We newly formulate the TBA equations for a defect with simultaneously non-vanishing 
reflection and transmission amplitudes. We indicate how these equations can be used to 
compute various thermodynamic quantities, which are, however, most interesting only when 
considered per unit length. By means of the TBA we compute the density distribution 
functions and use them to evaluate the Landauer conductance formula ([O]) for various defects 
in a complex free fermionic theory. We predict analytically the most prominent features in 
the conductance as a function of the temperature, i.e. the plateaux. 



We formulate the Kubo formula || for a situation in which defects are present (1.2) 



We evaluate the current-current correlation functions occurring in there by means of another 
non-perturbative method based on integrability, namely the bootstrap form factor approach 



[11, 12]. We provide closed formulae which solve explicitly the defect recursive equations 
involving any arbitrary number of particles. As for the Landauer formula, we also predict in 
this case the plateaux in the conductance as a function of the temperature analytically. 

There are several interesting open issues. Most challenging is to treat in full generality 
the massive and temperature dependent case of (|1.2|) . Unfortunately, the formulation of non- 
perturbative methods do not yet cover that situation [|5l] and it remains to be clarified how 
the form factor bootstrap program for the computation of two-point functions can be extended 
to that case. It would be further interesting to compute thermodynamic quantities per unit 
length by means of the TBA formulated in section 3.3. To classify possible defects more 
systematically is desirable even for free theories. 

Acknowledgments: We are grateful to the Deutsche Forschungsgemeinschaft (Sfb288), for 
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